Fractional exclusion statistics and the Random Matrix Boson Ensemble 
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The fc-body Gaussian Embedded Ensemble of Random Matrices is considered for N bosons dis- 
tributed on two single-particle levels. When k = N, the ensemble is equivalent to the Gaussian 
Orthogonal Ensemble (GOE), and when k — 2 it corresponds to the Two-body Random Ensemble 
(TBRE) for bosons. It is shown that the energy spectrum leads to a rank function which is of the 
form of a discrete generalized beta distribution. The same distribution is obtained assuming N 
non-interacting quasiparticles that obey the fractional exclusion statistics introduced by Haldane 
two decades ago. 
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I. INTRODUCTION 

Fractional exclusion statistics was introduced [l| as a 
generalization of the Pauli exclusion statistics. The par- 
ticles obeying such statistics are elementary excitations 
that can only exist in the interior of the whole system of 
interacting particles. It has also been shown that the ex- 
citations in exactly soluble many-body problems in one 
dimension, e.g. the Colagero- Sutherland model @, Q; 
the quasi-particles in the Luttinger model jj] and the 8 
function gas Q, obey this statistics. Excitations in the 
Haldane-Shastry spin model follow FES also £G| . In fact, 
it has rigorously been proved that quasi-particles with 
non-trivial exclusion statistics exist in a class of models 
that are solved by the Bcthe ansatz pj. Applications to 
other condensed matter systems have also been consid- 
ered 0. 

It has been shown [9| that a quantum liquid of par- 
ticles interacting via a long-range two-body potential, 
where the particles are supposed to be cither bosons or 
fermions, exhibits the nature of a quantum liquid with 
FES. The expectation that FES in arbitrary dimensions 
continuously interpolates between Fermi and Bose liquids 
holds only for one-dimensional systems. Recently [10(|, it 
has been shown that FES is a consequence of the inter- 
actions between the particles of the system and is due to 
the change from the description in terms of free-particle 
energies to the description in terms of quasiparticle en- 
ergies. 

On the other hand, several models in which two-level 
systems are occupied by bosons have been useful to un- 
derstand different quantum many-body systems and can 
be treated numerically by direct diagonalization. These 
models have as a common feature that the lower level has 
a scalar boson, but differ according to the multipolarity 
L of the second level. For example, the Lipkin-Meshkov- 
Glick model [ll[ in the Schwingcr representation has a 
scalar L = boson in the upper level, as is the case 
in the two-mode Bose-Hubbard model for Bose-Einstein 
condensates confined in two- well potentials [12( . A dipo- 
lar {L = 1) boson leads to the vibron model of quantum 



chemistry and a quadrupole (L = 2) boson corresponds 
to the interacting boson model [131 ] - 

Here we shall study a rather simple, yet general, two- 
level system with N spin-less bosons interacting through 
a /c-body random potential; this is called the fc-body 
Gaussian Embedded Ensemble for bosons [1J, [l5| . More 
specifically, we deal with the cases k = 2 and k = N. 
The former case corresponds to the most common physi- 
cal situation with two-body interactions. The latter cor- 
responds to the canonical ensembles of Random Matrix 
Theory (RMT) [16| . which is important in its own since 
RMT predicts accurately the universal statistical prop- 
erties of the fluctuations of the s pec tra of a large variety 
of complicated physical systems |17| ■ 

In this paper we use the fc-body boson random ensem- 
ble to calculate the integrated level density and show 
that, in all cases, beta-like rank distributions are ob- 
tained. Such beta-like distributions have been observed 
to hold in other complex systems in biology, social sci- 
ences and arts [18| . We also obtain here this distribution 
using quasiparticles that obey FES. This implies a gener- 
alization of two-level systems, since the second level must 
have a degeneration which is no longer an integer num- 
ber, as can occur with Haldane particles [l| as well as 
with Regge poles [19j . 



II. TWO-LEVEL BOSON SYSTEMS WITH 
RANDOM INTERACTIONS 



The fc-body Gaussian Embedded Ensemble for bosons 
is defined as follows [1J, [l5|, |2(|. We denote by k the 
rank of the interaction and by N the number of spin- 
less (scalar) bosons so k < N. The N bosons are dis- 
tributed in two single-particle levels; for simplicity, we 
assume that these levels are degenerate. We introduced 
the boson creation and annihilation operators a- and a, 
(j = 1,2), which satisfy the usual commutation rela- 
tions for bosons. The normalized A-boson states are 
denoted by |/^°) = (A i ,Jv)" 1 (a t 1 ) n (a^) Ar -"|0}, where 



A n ,N = [nl(N — n)!] 1 / 2 is a normalization constant and 
|0) is the vacuum state. Clearly, the dimension of the 
Hilbcrt space is Wb = N + 1 (sec Eq. (j4]) below with 
g = 0). Then, in second-quantized form, the most gen- 
eral Hamiltonian H^. involving fc-body interactions is 
given by 



Hf = 



E 

n.77i— 



Uk-r 



m (al)"(qT) fc -"(a 1 ) m (a 2 ) 
™' m A i. A i 



\k-r 



(i) 



r(/3) 



Stochasticity is built into the Hamiltonian -ffj, at the 
level of the fc-body matrix elements Vn,m- These matrix 
elements are assumed to be Gaussian distributed inde- 
pendent random variables, with zero mean and constant 
variance Vq = 1. As in the case of RMT [131, Dyson's pa- 
rameter j3 distinguishes the ensembles according to its 
symmetry properties: j3 = 1 corresponds to the case 
where time-reversal symmetry holds, and (3 = 2 holds 
when time-reversal invariancc is broken. The fc-body in- 
teraction matrix v^ ' is thus a member of the Gaussian 
Orthogonal Ensemble (GOE) for /3 = 1 or of the Gaus- 
sian Unitary Ensemble (GUE) for = 2 [l7|. 

We note that the number of independent matrix ele- 
ments of the fc-body matrix is equal to Kp(k) = /3(fc + 
l)(fc + 1 + 5/3,i)/2. Therefore, for fc < N, which is phys- 
ically the relevant case, the matrix elements of Hjf are 
correlated, i.e., the number of independent random vari- 
ables is much smaller than the number of independent 
matrix elements. In addition the matrix is sparse since 
many matrix elements are identically zero. The combi- 
natorial factors A n ^k and A m ,k are introduced in the def- 
inition of the ensemble Eq. ([T]) in order to have for k = N 
perfect identification, i.e., a one-to-one correspondence, 
with the canonical ensembles of RMT [1J, [15j • 

The numerical results presented below are related to 
the integrated level density, also referred to as the stair- 
case function, of individual realizations of the ensem- 
ble |T]). We take into account a large number of such 
realizations to achieve small statistical errors. For each 
member of the ensemble wc construct the staircase func- 
tion directly from the diagonalization of the Hamilto- 
nian. The staircase function is defined as Vl(E) = 

B(E - HJf } ) = Y,i@{E-Ei), where E, are the many- 
body eigenvalues of the specific realization of the Hamil- 
tonian; Q(x) is the Heaviside step function. By defini- 
tion, the staircase function is a monotonic non-decreasing 
function of the energy. For later convenience we shall use 
the rank function r(E) = TV + 1 — y\(E). Wc invert it 
numerically and express the energy E in terms of the 
rank r; the function E(r) is therefore a monotonic de- 
creasing function of the rank r. For each member of the 
ensemble, the rank distribution E(r) is fitted with the 
two-parameter distribution |l8[: 



f(r) = A 



(r n 



1-rf 



(2) 
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Here, A is the normalization constant of the distribution, 



FIG. 1: (a) Scatter plot of the fitting parameters a and b of 
Eq. ([2} computed for each realization of the interacting two- 
level bosons ensemble for k — 2 and N = 2000. The frequency 
distributions of a, b, and R 2 are given in figures (b), (c), and 
(d), respectively. 
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FIG. 2: Same as Fig. U for k = N. 



r niax = N + 1 is the maximum value attained by the rank 
r, and the constants a and b are two fitting exponents. 
We refer to f(r) as the discrete generalized beta distri- 
bution (DGBD). 

We present now numerical results on the fitting expo- 
nents a and b of Eq. [2j These were obtained by consid- 
ering 10000 members of the fc-body Gaussian Embedded 
Ensemble for two-level bosons, where we fixed TV = 2000 
bosons and /3 = 1; similar results were obtained for j3 = 2. 
The best fit was determined by a log-log multiple linear 
regression; the quality of the fit is measured by the square 
correlation coefficient R 2 , which must be close to 1 for 
a reliable fit. We focus on the cases fc = 2 and k = N 
since these are physically the most relevant. In Fig. Q] we 
present a scatter plot of the fitting parameters obtained 
for each individual member of the ensemble for k = 2 
and the frequency distributions of a, b, and R 2 ; Fig. [2] 
displays the corresponding results for k = N. Interest- 
ingly, these figures show that the distributions for a and 
b are much broader for k = 2 than for k = N; the results 
on R 2 simply give confidence to the fitting procedure. 

From the numerical results, we compute the average 
values of the fitting parameters. For k = 2 we have 
Sfc=2 = -6.43 x 10" 3 ± 0.1094 and b k=2 = 0.922 ± 0.209, 
where the assigned error corresponds to the standard 
deviation of the sample. Likewise, for k = N we ob- 
tain a k=N = -8.92 x 10~ 4 ± 2.15 x 10~ 3 and b k=N = 
0.924 ± 6.17 x 10~ 3 . In both cases, a is slightly nega- 
tive and very close to zero, and b < 1. We shall show in 
Section Hill that the fact that a ~ and b < 1 indicates 
that interactions mimic FES, i.e., interacting particles 
arc equivalent to non-interacting Haldanc quasi-particlcs 
with respect to the rank distribution. 



III. RANK DISTRIBUTIONS FOR 
FRACTIONAL EXCLUSION STATISTICS 

The combinatorial formula for the number of many- 
body states W g of N identical particles following FES 
and occupying a group of p states is given by (cf. page 
147 of [Hi) 



W a 



[p+(N-l)(l-g)}\ 
N\\p-gN-(l-g)}V 



(3) 



This expression reproduces for g = the well-known ex- 
pression for bosons, Wr = [p + N — 1]!/[AM (p — 1)!], and 
for .9 = 1 that of fermions W F = p\/[N\ (p - N)\}. Note 
that Eq. ([3]) can be rewritten as 



W g = 



\p' + JV-1]! 
N\[p' - 1]! 



with 



p'=p-g(N-l). 



(4) 



(5) 



One has therefore N identical bosons occupying an ef- 
fective group of p' single-particle states. Note that p' 



could be formally negative. This gives the possibility of 
an interesting connection with the famous Regge trajec- 
tories used in high-energy physics [211 ] . 

In order to compare the rank distribution of iV par- 
ticles obeying FES with those obtained in Section |TT] 
for bosons, we shall group the p' single-particle levels 
filled by N Haldane particles into two levels with ener- 
gies eo = and ei > eg, with p' = p' + p[. If p' Q = 1, 
the number of (non-interacting Haldane particles) mi- 
crostates with energy E = n\€x is given by 



w(nl = E/ei) 



(m+pi-i)! 



(6) 



m!(pi-i)! 

The rank distribution with respect to the energy is then 



r(E) 
which for p' = 1 yields 



N 

E 



u(v), 



(7) 



r(E) 



N 

E 
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(iV + Pi)! QE-ei+pi)! 
iV! Pl ! (£-ei-l)!pi! 



(8) 



Taking the continuum limit of E and using Stirling for- 
mula, we obtain 



E~(N P '^ -p'^.r) 1 '^ 



(9) 



For particles obeying FES the effective number of 
states Pi can be positive or negative, as we have remarked 
in connection with Eq. ([5]) . If p[ > and N is large, we 
have 



E 



— ^ = (7V-r) 1/p 'i 



(10) 



with A/" = N Pl /pi\, which is a discrete generalized beta 
distribution with a = and b = l/p^, as seen from direct 
comparison with f(r) in Eq. ^). On the other hand, if 
p[ < and N is large, we have a rank distribution of the 
form 



[(-Pi) 



iVrf-i/Pi 



(11) 



which corresponds to the values a = l/\p[\ and b = of 
the DGBD. In order to have (—p[)\ real, p' x must satisfy 
— 2,s — 1 < pi < —2s, with s a positive integer. 



IV. DISCUSSION 



We have given numerical evidence to show that the 
fc-body embedded random matrix ensemble for bosons 
leads to an integrated level density which corresponds to 
a discrete generalized beta distribution. We have defined 
the ensemble for N bosons interacting through a fc-body 
force distributed in two single-particle levels. When k = 
N the GOE is obtained by construction, and the two- 
body random ensemble 22] for bosons corresponds to 
k = 2, a more physical case. 

The discrete generalized beta distribution fits rather 
well the statistical regularities in many instances: in par- 
ticular, for the rank-citation profile of scientists [23| . This 
distribution is present in systems whose output depends 
on the restricted difference of statistical variables j24| . 
Examples of applicability are neuron dynamics (excita- 
tion and inhibition) and ecological community evolution 
(birth and death). 

We then considered a set of independent Haldane 
quasiparticles, which obey fractional exclusion statistics, 
distributed in the same two-level system. It is shown that 
DGBD is also obtained. With respect to this distribu- 
tion, at least, the above results indicate that the interac- 
tion among the particles makes the system equivalent to 
a set of quasiparticles obeying FES. 
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